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The leading order contribution to the shear and bulk viscosities, ry and i^, of a gluon plasma 

in perturbative QCD includes the gg -f-)- gg (22) process, gg -f-)- ggg (23) process and multiple 

^SJ ■ scattering processes known as the Landau-Pomeranchuk-Migdal (LPM) effect. Complete leading 

T-H ■ order computations for rj and ^ were obtained by Arnold, Moore and Yaffe (AMY) and Arnold, 

^^ ■ Dogan and Moore (ADM), respectively, with the inelastic processes computed by an effective g -^ gg 

O^ , gluon splitting. We study how complementary calculations with 22 and 23 processes and a simple 

^ ■ treatment to model the LPM effect compare with the results of AMY and ADM. We find that 

r- { ' we agree with their results within errors. By studying the 23 contribution to r), we find that the 

^^ , minimum angle among the final state gluons 9 in the fluid local rest frame has a distribution that is 

peaked at ~ y/ci2, analogous to the near coUinear splitting asserted by AMY and ADM. However, 

the average of 9 is much bigger than its peak value, as its distribution is skewed with a long tail. 

The same 9 behavior is also seen if the 23 matrix element is taken to the soft gluon bremsstrahlung 

limit in the center-of-mass (CM) frame. This suggests that the soft gluon bremsstrahlung in the 

f~ l ■ CM frame still has some near coUinear behavior in the fluid local rest frame. We also generalize our 

^ I result to a general SU{Nc) pure gauge theory and summarize the current viscosity computations in 

' QCD. 

D. 

I \. I. INTRODUCTION 

^ ■ Shear and bulk viscosities, rj and C, are transport coefRcients characterizing how fast a system goes back to equi- 

fSJ ', libriuni under a shear mode perturbation and a uniform expansion, respectively. In a weakly interacting hot gluon 

CnI ■ plasma, rj is inversely proportional to the scattering rate, 770c l/F ex 1/a^ InaJ^ [ll, where as is the strong coupling 

constant. ( is suppressed by an additional factor of (TT) , arising from the response of the trace of the energy 

momentum tensor [T^j to a uniform expansion. Thus, ^ vanishes when the system is "conformal" or scale invariant. 

For a gluon plasma, the running of the coupling constant breaks the scale invariance. Thus, T*^ oc /?(as) c>c a^, 

(^ oc a'^/lnaj^ [2]. In the perturbative region, C/77 ex a^ ^ 1. 

T— I ■ In the strong coupling region, smaller 77 is expected. The so-called "perfect fluid" is a fluid with the smallest shear 

^^ I viscosity per entropy density (s) ratio, 77/s. It is conjectured that rj/s has a minimum bound l/(47r) [3j. This is 

. ^H • motivated by the uncertainty principle of quantum mechanics because 77/s is related to AEAt, the mean energy and 

^^ ' life time of quasiparticles. While the number l/(47r) arises from the universal value rj/s = l/(47r) obtained for a big 

H . class of strongly interacting conformal field theories (CFT's) in the large N, N being the size of the gauge group, 

- - ■ and large t'Hooft coupling limits [3-[a|- This class of strongly interacting CFT's are dual to another class of weakly 

interacting gravitational theories in anti-de-Sitter space backgrounds. This anti-de-Sitter space/conformal field theory 

correspondence (AdS/CFT) (a-Q allows that 77/s in strongly interacting CFT's can be computed in weakly interacting 

gravitational theories. 

The smallest 77/s known so far is realized in the hot and dense matter (thought to be a quark gluon plasma of 
QCD) just above the phase transition temperature (Tc) produced at RHIC |9l4ll| with rj/s = 0.1 ± O.l(thcory) ± 
0.08(experiment) [l^l- A robust upper limit 77/s < 5 x l/(47r) was extracted by another group [13 and a lattice 
computation of gluon plasma yields 7;/s — 0.102(56) (at temperature T = 1.24Tc) ^^. Away from Tc, rj/s of QCD 
becomes larger due to small couplings at high T or small derivative goldstone boson couplings at low T. We will 
summarize the current status of QCD 77/s vs. T in Fig. [71 

As for the bulk viscosity C,, it is small in the perturbative region. However, near T^ the rapid change of degrees of 
freedom in T gives a rapid change of T'' which could give very large C/s [l5l.[l6|. 

The best perturbative QCD calculation of C was carried out by Arnold, Dogan and Moore (ADM) [2| using the 
same approach as the -q computed by Arnold, Moore and Yaffe (AMY) in Refs. [ll, ll3|- In both 77 and C, the leading 
order (LO) contribution involves the elastic gg <-^ gg (22), inelastic number changing process gg O ggg (23) and 
multiple scattering processes known as the Landau-Pomeranchuk-Migdal (LPM) effect. In the complete leading order 
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computations for 77 and C obtained by AMY and ADM, respectively, the inelastic processes were computed using an 
effective g ^ gg gluon splitting obtained after solving sophisticated integral equations. 

In this paper, we study how simpler calculations with the 22 and 23 processes and a simple treatment to model 
the LPM effect compare with the results of AMY and ADM. This approach is similar to the one used by Xu and 
Greiner (XG) |18l[l9| who claimed that the dominant contribution to 77 is 23 instead of 22, in sharp contradiction to 
the result of AMY. If XG's assertion is correct, then the small value of rj/s measured by RHIC can be explained by 
the conventional perturbative QCD picture alone. 

While our approach is not model independent because our simplified treatment of the LPM effect, it can be used 
to double check XG's result since the two approaches are very similar. We find that we cannot reproduce XG's result 
unless the 23 collision rate is at least multiplied by a factor 6 (part of this result was asserted in Ref. [20|)- In the 
mean time, our 77 agrees with AMY's within errors, while our C also agrees with ADM's within errors. 

Although our result does not provide a model independent check to AMY and ADM's results, we can still study 
the angular correlation between final state gluons using our approach. Because the 23 matrix element that we use is 
exact in vacuum, we can check, modulo some model dependent medium effect, whether the correlation is dominated 
by the near collinear splitting as asserted by AMY and ADM. 

We study the distribution of the minimum angle 9 among the final state gluons. If the near collinear splittings 
dominate, then most probable configurations would be that two gluons whose angles are strongly correlated and their 
relative angle tends to be the smallest among the three relative angles in the final state. This can be seen most easily 
in the center-of-mass (CM) frame of the 23 collision with two gluons going along about the same direction while the 
third one is moving in the opposite direction. We expect it is also the case in the fluid local rest frame. 

We find that 6 is peaked at 6* ~ ^/oil, analogous to the near collinear splitting asserted by AMY and ADM. However, 
the average of 6, (0), is much bigger than its peak value, as its distribution is skewed with a long tail. 

The same 9 behavior is also seen if the 23 matrix element is taken to the soft gluon bremsstrahlung limit in the 
CM frame. This suggests that the soft gluon bremsstrahlung in the CM frame still has some near collinear behavior 
in the fluid local rest frame. 

We also generalize our result to a general SU{Nc) pure gauge theory and summarize the current viscosity compu- 
tations in QCD. 

II. KINETIC THEORY WITH THE 22 AND 23 PROCESSES 

In this section, we will focus on the ( computation. We refer the formulation for calculating 77 to Ref. f2(j|. 
Using the Kubo formula, C can be calculated through the linearized response function of a thermal equilibrium 
state \il) 

C = hnr ^ r dt I Sr e^^' {n I [T^ (x) , T^ (0)] \n) . (I) 



>o 9a; 7o 

In the LO expansion of the coupling constant, the computation involves an infinite number of diagrams [2ll, |22| |. 
However, it is proven that the summation of the LO diagrams in a weakly coupled (p"^ theory ^21-25] or in hot 
QED [26i| is equivalent to solving the linearized Boltzmann equation with temperature-dependent particle masses and 
scattering amplitudes. This conclusion is expected to hold in perturbative QCD as well. 

The Boltzmann equation of a hot gluon plasma describes the evolution of the color and spin averaged gluon 
distribution function fp(x) which is a function of space-time x = (i, x) and momentum p = {Ep, p). 

The Boltzmann equation for the gluon plasma |27h32| | reads 
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^v^"^" " ^.S^(-'0yi..,„-/'^^ 



(2) 



The collision kernel 
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c?ri...,^(i+i)...(„_i)p = ]J .^^.3^^^^ ^|Mi...;^(;+i)...(„_i)p|^(27r)^,5^(^p^ - ^ v,-v) (3) 

has summed over all colors and helicities of the initial and final states in the matrix element squared. Ng = 2{N^ ~ 1) = 
16 is the color {Nc = 3) and helicity degeneracy of a gluon. The i-th gluon is labeled as i while the 77-th gluon is 



labeled as p. For a process with / initial and (n — /) final gluons, the symmetry factor N{n, I) = l\{n — I — 1)!. For 
example, processes 12 — >■ 3p, 12 — > 34p, 123 — ?► Ap yield (n, ?) = (4,2), (5,2), (5,3) and N{n,l) = 2,4,6, respectively. 
|A/i...;_>(/4.i)...(„_i)pP is the matrix element squared for the process 1 • • • Z — )■ (Z + 1) • • • (n — l)p without average over 
the degrees of freedom for incident gluons, i.e. it includes a factor N'^. 



In vacuum, the matrix element squared for the 22 process is 



|M: 



12->34| 



SNgiAira.N, 



tu su st 

g2 f2 y2 



(4) 



where as — g'^/{4:Tr) is the strong coupling constant, and {s,t,u) are the Mandelstam variables s = {pi +^2)^, 
t = (Pi -Pa)^ and u = {pi - Pi^ ■ 



my 



For the 23 process J33l-l35l|. under the convention X)i=i-P« — Oj "^^ have 



|Afl2345-i-o| — |A'^0^12345| 
1 



-iV,(4™,7V,)3 [(12)^ + (13)^ + (14)^ + (15)^ 

+ (24)* + (25)* + (34)* + (35)* + (45)*^ 

1 



E 



porm{l,2,3,4,5} 



(12) (23) (34) (45) (51) 



(5) 



where {ij) = pi ■ pj and the sum is over all permutations of {1,2,3,4,5}. To convert to the convention pi + P2 — 
P3 + Pi + P5, we just perform the replacement: 



|-Wl2->345| — |Afo^l2345| 



\M: 



345^121 



\M: 



12345^01 



Pi^-pi,P2->-P2 



Pl->-pl,P2->--P2 



(6) 



In the medium, the gluon thermal mass serves as the infrared (IR) cut-off to regularize IR sensitive observables. 
The most singular part of Eq.((31) comes from the collinear region (i.e. either t sa or u w 0) which can be regularized 
by the HTL corrections to the gluon propagators |36l . |37J and yields [S^ , 



|Mi2^34|' « AiATrasN,fNg{AEiE2)' 



(1 — x'^) COS( 



Hl q2(l-a;^ 



n7 



(7) 



where q = P2 — Pa — {qo, q), x — qo/\q\ and 4> is the angle between Pi x q and p2 x q. The HTL self-energies 11^ 
(longitudinal) and IIt (transverse) are given by 



IIl = m^ 



n7 



X 1 + X n_ 

1 In h i—x 

2 1-x 2 



h -(1 - X ) In i-xil - x^) 

2 V ' l-x 4 ^ ' 



(8) 



The external gluon mass rrioo (i-e. the asymptotic mass) is the mass for an on-shell transverse gluon. In both the 
HTL approximation and the full one-loop result, m^ — H^ {\x\ = 1) = jrijj/2, where niD — {^nasNc/SY^^T is the 
Debye mass. 

Previous perturbative analyses showed that the most important plasma effects are the thermal masses '^ gT acquired 
by the hard thermal particles |39l - l4l| . So a simpler (though less accurate) treatment for the regulator is to insert inn 
to the gluon propagator such that 



|Afi2^34|'«8iV<,(47ra.iVcs)' 



1 



1 



(i - m]jY [u ~ ml) 



2 \2 



(9) 



It can be shown easily that Eqs. Q and ([9|) coincide in the center-of-mass (CM) frame in vacuum. This treatment 
was used in Refs. [Mliiliil. 

Eq. ([9]) is often expressed in qr, the transverse component of q with respect to pi, in the CM frame. If we just 
include the final state phase space of the i-channel, near forward angle scatterings (q^ w q^ w 0), then the backward 
angle contribution from the u-channel can be included by multiplying the prefactor by a factor 2 

o2 



\Mi2^34\cM 



16Ng{4nasN,y 



(q^+m^) 



2 \2 ■ 



(10) 



But if one includes the whole phase space in the calculation, then the factor 2 is not needed: 

\M,2^,,\Im ~ ^N,{A^asN,f ' ^ (11) 

Note that the constraint q^ w is removed because both the near forward and backward scatterings have small c^ 
but only the near forward scatterings have small q^. 

For the 23 process, because the matrix element is already quite complicated, we will just take ttid as the internal 
gluon mass as was done in the 77 computation in Ref. [20'] and then estimate the errors. In the X]i=i Pi = convention, 
one can easily show that an internal gluon will have a momentum of ±(pi +Pj) rather than ±(pi ~ Pj)- Therefore, 
the gluon propagator factors (ij) in the denominator of Eq. ([5]), should be modified to 



Pt ■ P] 



2m^ - m\, 



= Pi -P] , (12) 

where we use the asymptotic mass for the external gluon mass. Then one applies Eq. ^ for the Boltzmann equation. 
In the numerator, the {ijY combination is set by T and is 0{T^). So we can still apply the substitution of Ea. (|12l) . 
even if the {ij) factors might not have the inverse propagator form. The error is ^ Tnjjiij)^ — 0{asT^), which is 
higher order in as- 

It is instructive to show that Eqs. (|5l6p and p^ give the correct soft bremsstrahlung limit. Using the light-cone 
variable 

P = {p^,P^,Pt) 

= {pq+P3,Pq-P3,Pi,P2) , (13) 

we can rewrite one momentum configuration in the CM frame in terms of p,p', q and k: pi ^ p, p2 = p' , ps = p + q — fc. 
Pi — p' — q and p^ — k, with 



The on-shell condition p§ — pi 





P = (V7,TO^/V^, 


0,0), 






P = ("tL/%/S: \/s. 


0,0), 






k = {yVs, {kr + ml^) /2/Vs, 


kT,0), 




q = {q^,q^,(iT)- 






ml^ 


yields 






?+ 


~ -g|./^/i, 






q" 


fc|, + yq^ - 2ykT 


Qt + (1 - 


V + V^)ml^ 



(14) 



2/(i-y) 



(15) 



Here y ~ k'^ jp'^ — kTC^ I \fs is the light-cone momentum fraction of the bremsstrahlung gluon and z is its rapidity. 
In the central rapidity for the bremsstrahlung gluon, i.e. 2; ~ 0, y is toward zero. In this case ps = fc is very small 
compared to p\ and pi . 

Now, in the limit s — ?> cx), y — ^ 0, while keeping y^/s fixed, we have 

p = (v^,o,o,o), 

p' = (0,^,0,0), 
k = (yVs, (fc|-Hm^)/2/v^,fcT,0), 
q = (0,(4-FTOL)/yV^,qT)- (16) 

In this limit, pi-4 are hard (their three momenta are O (\/s)) while q{= P2 — Pi = —pi + P3 + P5) and p^ — k are 
soft (their three momenta are much smaller than ^/s) . In this particular limit of the phase space, the matrix element 
becomes 



|Mi2->345|cM ~ ^HAnasNafNg 



s^ 



{k^ + ml^) (q^ + 171%) {kT-qr^+ml, 



(17) 



where the prefactor is equivalent to 34567r'^Q;g7V^ when iVc = 3. 6 different permutations of (^3,^4,^5) give the 
same expression as Eq. (J17p due to the permutation symmetry of Eq. ([5]) . Those permutations are corresponding to 
different symmetric diagrams, just as the two permutations of (^3,^4) in Eq. ([9]) give the t- and u-channel diagrams by 
the crossing symmetry. Analogous to Eqs. pUj) and ([TT|l . if we only include the constraint phase space of (^3,^4,^5), 
then we need to multiple Eq. p7p by a factor 6 to take into account the permutations of (^3,^4,^5). But if we include 
all the phase space in the calculation, then Eq. (fT7|) will be the one to use. Any additional symmetry factor will result 
in multiple counting. 

The ratio of Eq. ^T7\ to Eq. PT|) reproduces the Gunion-Bertsch (GB) formula [45] after taking tod,?71oo -> 0. 
One can find the derivation of the GB formula from Eq. ([SJ in Appendi x Rl One can also expand the "exact" matrix 
element in Eqs. (|5l6p in terms oit/s to extend the GB formula f46ll48l. l4g| . 

An intuitive explanation of the LPM effect was given in Ref . [43] : for the soft bremsstrahlung gluon with transverse 
momentum kx, the mother gluon has a transverse momentum uncertainty ^ fc^ and a size uncertainty ~ l/fcr- It 
takes the bremsstrahlung gluon the formation time i ~ 1/ (fcr^T) ~ -E'/t/fcr^ to fly far enough from the mother gluon 
to be resolved as a radiation. But if the formation time is longer than the mean free path Imfp ~ 0{a^^), then the 
radiation is incomplete and it would be resolved as gg — >■ gg instead of gg — > ggg. Thus, the resolution scale is set by 
t < Imfp- This yields an IR cut-off fc|, > Ek/lmfp ~ 0{as) on the phase space \44^. Thus, the LPM effect reduces the 
23 collision rate and will increase 77 and C. Our previous calculation on rj using the Gunion-Bertsch formula shows 
that implementing the ttid regulator gives a very close result to the LPM effect [i^]- Thus, we will estimate the size 
of the LPM effect by increasing the external gluon mass rUg from moo to mo- 

III. AN ALGORITHM BEYOND VARIATION TO SOLVE FOR C 

Following the derivation of Ref. |22l ] , the energy momentum tensor of the weakly interacting gluon plasma in kinetic 
theory can be modified as 



d^p 



T^,uix) = Ng I 7:^:T3Tr/p(a;) (PpP^^ - U (x) g^,^) , (18) 



where the U {x) term is an effective mass term from the medium effect and Ep = -^/p^ + rn^ . When the system 
deviates from thermal equilibrium infinitesimally, fp{x) deviates from its equilibrium value f!:^ = (e"'^'^ — 1)^^ 



fp = !;"• + 5 fp. (19) 



And so does Tn^: 



-^^A- = ^9 / TTT^SU P,Pu - t'M^.T^— |i , (20) 



where the energy momentum conservation d^T^^i, — has been imposed. 

In hydrodynamics, small deviations from thermal equilibrium can be systematically described by derivative expan- 
sions of hydrodynamical variables with respect to spacetime. We will be working at the v(a;) — frame for a specific 
spacetime point x (i.e. the local fluid rest frame). This implies d^v^ = after taking a derivative on v^{x)v^^{x) = 1. 
Then energy momentum conservation and thermal dynamic relations (we have used the property that there is no 
conserved charge in the system) in equilibrium allow us to express the time derivatives dtT and dtv in terms of the 
spacial derivatives V • v and VT. Thus, to the first derivative expansion of the hydrodynamical variables v and T, 
the bulk and shear viscosities are defined by the small deviation of T^^i, away from equilibrium: 



dv' dvi 2 
dx^ dx^ 3 



5T„- = -C%V . V - ry — + — - -L,V . V ) , (21) 



where i and j are spacial indexes. Also, 5T{)i{x) = 0, since the momentum density at point x is zero in the local fluid 
rest frame, and one defines Too to be the energy density in this frame. Therefore, 

JToo = ^NgJ j0^Sfp (p2 + m2) , (22) 



where 



fh'^ml-T'^ = -iiV,/3(g2)T2 = ^^N^a^T^. (23) 



Matching kinetic theory (Ea. ipO]) ') to hydrodynamics (Eq. ipT]) ) to the first derivative order, 5fp can be parametrized 



as 



where 






su = -xp.f;' {I + /;'), 

B,j{p) I ( dv' ^ dv' 2 



T 2 \dx^ ' dx^ 3%^-^ 



(24) 



(25) 



We can further write Bij{p) — B{p){p^p^ — ^Sij) with p the unit vector in the p direction. A{p) and B{p) are 
functions of p. They can be determined by the Boltzmann equation to give the solution of the bulk and shear 
viscosities, respectively. In this work, we will focus on solving the bulk viscosity. 

Working to the first derivative order, the Boltzmann equation becomes a linear equation in Sfp which yields 

^ - clip' + m') - 2^ / rfri2^3p/rV2'(l + /r)(/r)-M^3 + Ap-A,- A,] 

Q J 123 

+^ I rfri2^34p/r/2'(i + /r)(i + fi')u;')-'[A-i +a,+Ap-a,- a^] 

^^^g J 1234 

+^ / dri23^4p/rvr/r(i + /r)(/;'')-M^4 +Ap-a,^a^- as]. (26) 

D-'^S J1234 

Here we have used the notation Ap = A(p) and Ai = A{pi) with i = 1, 2, 3. The speed of sound squared cj. is defined 
as [3,[2l, 



dP _ /rfVp"^(l + /p^^)p' 



de i!d^pfp\i + f;'){p' + . 



Then Eqs. ()20I21I25|) yield 



^ T J {2ttYEp^p y ^Jp ) 



3P'-c?(p'+™') 



Aip), 



(27) 



(28) 



where we have added the c^ term which is proportional to ST^q = for convenience. By substituting Eq. (j26l) into 
Eq. dlS]), we obtain 



c-5^/n 



(Pki 



8T J ^^ (2tt)32E, 



|Mi2^34p(2^)^<54(£;i +E2-Es-Ei) 



X (1 + ./r)(l + r2')ft''fr[A^ + A,-A,- A^f 



19T / 11 



d^k, 



12T J -•■-•- (27r)32K 

F'^9\/^ 1 feq\feq feq feq 



Mi2^345r(27r)*<5*(^i +E2- E3 -Ei-E5 



X (1 + /r«)(l + fl')fl'frr,''[A, +A, + A,-A,- A2Y. 



By the definition of Cg, the following integral vanishes: 

(27r)3£;p-'P ^ ■'P 



cj(p'+m^) 



Sp = 0. 



(29) 



(30) 



We will use this property later. 

Now we first review the arguments that cast the computation of C as a variational problem [171 . |50| . Then we show 
how one can go beyond variation to find the answer systematically. Let us rewrite Eq. (j26p schematically as 



and Eqs. (|28l29p as 



\S)^C\A), 
C^{A\S) = {A\C\A) 



(31) 
(32) 



Note that Eq. ^ is just a projection of Eq. §^. Using \A) = C'^ \S) from Eq. dHD, 

C=(5'|C-i|S'). (33) 

Technically, finding an ansatz Aam that solves the projected equation {S\Aanz) = {Aanz \C\ Aanz) of ([32]) is easier 
than solving the original integral equation (PT|) . But this will not give the correct viscosity if C \Aanz) 7^ \S). However, 
the resulting bulk viscosity has the property 

Qanz = ^ \Aanz \^\ Aanz) + ^ \Aanz\^) 
= ^ \^anz\^\ ^anz) ~^ \S \^ \ S) 

< {S\C-'\S)^C, (34) 

where l^^n^) = \Aanz) — C^^ \S) and (^Jm^ \C\ A'^nz) is real and non-negative. Thus, a variational calculation of C, is 
possible: one just demands {S\Aanz) = {Aanz \C\ Aanz) and try to find the maximum Qanz- In what follows, we show 
an algorithm that will approach the maximum Qanz systematically. 

We will choose a basis {Ai\i = 1, 2, . . . , n} with n orthonormal functions satisfying 

{A..\C\A,)^5,,. (35) 

Also, we impose 

j0^/p'(^ + /P')(P' + f^')Mp) = 0, (36) 

such that when we take the ansatz for A 



4nUE^^^- (37) 



1=1 



the constraint STqq = is automatically satisfied. Then Eq. ^32\i yields 



n 



ciil = E ^^ (^^1^) = E '^^^^^ (^^ ici ^.) = E ^- (38) 

2—1 ^J — 1 "i— 1 



The above equation can be rewritten as 

n 



n 
\2 



= E(^»i^)'-E(^''-(^^i^))'- (39) 



i=l i=l 



Thus, the solution 



rf. = (A, 1 5) (40) 

satisfies the projected equation f|38p and maximizes d^anz as well. So we have 

n 

dl-E(^^i^)'- (41) 

i=l 

Then since (viols') is real, Canz is monotonically increasing with respect to n. 
Below are some comments about this algorithm: 

(I) Canz < Canz < Canz °° = C- This algorithm Systematically approaches C, from below. 

(II) This approach solves for a vector Aanz satisfying (A[mz\S) = (^inz |C'|v4)i„^ V The solution is independent 



of the basis used. One can show that all bases which are linear combinations of Ai give the same Cani = ( S\A. 



We will use the following basis: 

n 

A(")(p)=^c,(i?p/Tr. (42) 

i=0 

We have used the property (II) discussed above such that we need not make the basis orthonormal. The Ep term will 
not contribute to (A'^") |C|^^"^) or <^S'|^("^) because of energy conservation in Eq. (P^ and because of Eq. ([50)) . 
But it can be used to make the constraint STqo = satisfied. 
An alternative basis is used in Ref. j2|: 

^''-"'iP^-t/^ ^jT+lr-^ +d'E,- (43) 

The two bases give give consistent C. For example, at a^ = 0.1, the agreement is better than 1% when we work up to 
n = 6. 

IV. Nc SCALING AND NUMERICAL RESULTS 

A. Nc Scaling 

Viscosities of a general SU{Nc) pure gauge theory can be obtained by simply rescaling the SU{3) result. Using the 
above formulas, it is easy to show that 

C = Ngg, (a,7Ve) T\ r, = N^g^ {a,N,) T\ (44) 

where gi and (72 are dimensionless functions of agNc only. This, together with s oc Ng, yields 

^ = /ii (a.iVe) , - = /i2 (a.iVe) , (45) 

s s 

where hi and /i2 are also dimensionless functions of UgNc only. Thus, our C/s, 77/s and Q/rj vs. agNc curves in Fig. [6] 
are universal and suitable for a general SU{Nc) pure gauge theory. From now on, Nc ^ 3 unless otherwise specified. 
One can always rescale the results to an arbitrary Nc- 

B. Leading-Log result 

As discussed above, in the leading-log approximation, one just needs to focus on the small qt contribution from 
the 22 process while setting co = 0. Furthermore, it was shown in (Sll, l52| that using the HTL regulator ([7|) gives the 
same LL shear viscosity to that using the mn regulator ^. For the bulk viscosity, this is also true. We obtained the 
same LL result as |2|, 



^■'■■'"■''HWr (*' 



This can be compared with [l|, |20| 



ai\n{l/g)- 



For a gluon plasma, we have 



'?"^0-17^:j177T77T- (47) 

^~2.6a^ = 48(l/3-c2)'. (48) 

VLL ^ ' 

This is parametrically the same as (^/rj = 15 (1/3 — c^) for the absorption and emission of light quanta (e.g. photons, 
gravitons or neutrinos) by the medium [53]. In the as <S^ I region where QCD is perturbative, C ^ V- Using the 
entropy density for non-interacting gluons, s — Ng'^jr-T^ , we have 

^.0.063-^, ^-4^. (49) 

s ln(l/g) s al\TL{l/g) 



Figure 1: Left panel: The ratio of our numerical result (denoted as "exact") for the shear viscosity i] to AMY's. The error 
is bounded by the upper and lower bound. The result using the GB matrix element [Eq. (|17|l ] is also shown. Right panel: 
Comparison of our result of rj/ri22 with those of AMY, XG and "GBx6"(see the text). For our result the HTL gluon propagator 
is used for the 22 process, the "exact" matrix element Eqs. (|5I6I I12p are used for the 23 process, and the external gluon mass 
is set to rrioo. In both panels, the error bands are shown. AMY's result is taken from Ref. [l|, [ij and XG's result is from Ref. 
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Figure 2: rj/s with various inputs. 
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C. Numerical results of rj and ^ 



111 our calculation, we use the HTL propagator for the 22 process. For the 23 process, f or technic al reasons, we 

~ \/p^ + m^ in kinematics 



use the internal gluon mass mo instead of the HTL propagator in Eqs. 
and /p^ for the external gluon distribution. The errors from not implementing HTL propagator in the 23 process 
and the modeling of the Landau-Pomeranchuk-Migdal (LPM) effect, and from the uncalculated 0(^07) higher order 
corrections are estimated in Appendix [Bl 

In Fig. [U we show our main result for the shear viscosity 77 in our previous paper [42] , together with the theoretical 
error band bounded by "upper bound" and "lower bound" curves. Note that previously we estimated the higher order 
effect to be 0{as) suppressed. But since the expansion parameter in finite temperature field theory is g instead of 
g^, we enlarge the error of the higher order effect to 0{y^) here. The result agrees with AMY's result within errors 
in the left panel although our central value is lower at larger as- If we replace the "exact" matrix element of Eqs. 
(I51 l7ll2p by the GB matrix element of Eq. ([T7]) , then ij is reduced by still close to the estimated lower bound. This 
means the 23 collision rate in GB is bigger than that in "exact." 

The effect of the 23 process can be seen more clearly in the ratio 77/7722 (?722 means the shear viscosity with the 
22 process included only) shown in the right panel, where we also show AMY's and XG's results for comparison. In 
AMY's result [ll, [ITI, the near collinear 1 -H> 2 process gives 77/7722 close to unity. This implies their 12 collision is just 
a small perturbation to the 22 rate. However, XG employ the soft gluon bremsstrahlung approximation in the matrix 
element for the 23 process, gives 7^/7/22 — 0.11 ~ 0.16 around 1/8 in Ref. [13, indicating that their 23 collision rate 
is about 7 times the 22 one. In their improved treatment using the Kubo relation ^19] , they give 77/7722 — 0.1 ~ 0.3, 
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Figure 3: Left panel: comparison of our result for the bulk viscosity (" and its error band (see the appendix) with ADM's result. 
'LL' denotes the leading-log result of EQ. (|46p . 'ADM LO' denotes ADM's leading order result read off from Ref. [2] (only 
available for Qs > 8 x 10"''). Right panel: our full result for (^ (denoted as 22+23) and ^ with the 23 process only (denoted as 
23). 
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indicating that the 23 collision rate is about 2~9 times the 22 rate. 

Our central result lies between AMY's and XG's results. However, even consider the lower bound, our 23 rate does 
not get bigger than the 22 rate. Thus, it is qualitatively consistent with AMY's result but inconsistent with XG's 
result. When compared with AMY's result, in addition to the error band shown in the left panel, there is still ^ 10% 
difference at a^ = 0.01. This is consistent with the ^ m'i/T'^ effect from using different inputs to external gluon 
mass — we use rUg while AMY uses zero. 

We find that we cannot reproduce XG's result unless we use a 23 matrix element squared at least 6 times larger. 
To compare with XG's calculation, we use the same ttIoo = and LPM effect as XG, and 



|Mi2^345lcA/-^6x54567V^ 



2 ' 






fc2, [q'^ + to|,) 


(kr- 


- qrf + m?jj 



(50) 



which is a slightly different variation of the GB matrix element squared of Eq. ([T7)) multiplied by a factor 6 (denoted 
as "GBx6"). This reproduces XG's result at a^ = 0.01. The origin of this discrepancy is yet to be resolved. 

In Fig. [21 ri/s with various inputs are shown. At as — 0.3 and 0.6, the GBx6 curve yields -q/s — 0.19 and 0.09, 
respectively, while XG has 0.13 and 0.08. The central value of the "exact" result is about two times lager. 

Our result for the bulk viscosity Q using the "exact" matrix element for the 23 process is shown in Fig. [3] We have 
worked up to n = 6 and seeir good convergence. For example, we obtain Q„2 ^ 95%, 98%, 99.5% of Q.nz at Ug — 10~^ 
for n — 3,4,5 respectively. The convergence for larger ag is eveir better. When a^ < 10~®, our result approaches 
the LL one. At larger a^, the 23 process becomes more important such that when a^ ^0.1, ^ is saturated by the 23 
contribution (see the right panel of Fig. [3]). Our result agrees with that of ADM [2| in the full range of as within the 
error band explained in the Appendix |B] 



D. Angular Correlation in 23 Process 



As mentioned in the introductioir, although our result does not provide a model independeirt check to AMY and 
ADM's results, we can still study the angular correlation between final state gluons using our approach. Because 
the 23 matrix element that we use is exact in vacuum, we can check, modulo some model dependent medium effect, 
whether the correlation is dominated by the near coUinear splittings as asserted by AMY and ADM. 

We study the distribution of the minimum angle 9 among the final state gluons. If the near collinear splittings 
dominate, then the most probable configuration would be that two gluons whose angles are strongly correlated and 
their relative angle tends to be the smallest among the three relative angles in the final state. This can be seen most 
easily in the center-of-mass (CM) frame of the 23 collision with two gluons going along about the same direction while 
the third one is moving in the opposite direction. We expect it is also the case iir the fluid local rest frame. 



We find that 9 is peaked at 6*^' 



eak 



fal, analogous to the near collinear splitting asserted by AMY and ADM. 



However, the average of 9, (9), is much bigger than its peak value, as its distribution is skewed with a long tail. Below 
are more detailed descriptions of our results. 
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Figure 4: The distribution of 6 in the fluid local rest frame. Left panel: weighted by the phase space and the Bose-Einstein 
distribution functions only. Middle panel: weighted by the contribution to r;23 with the "exact" matrix element. Right panel: 
weighted by the contribution to r;23 with the GB matrix element. 
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Figure 5: Same as Fig. |4]but with Q in the CM frame. 
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We show the distribution of B in the fluid local rest frame in Fig. SI and show the distribution in the CM frame 
of the 23 collision in Fig. [S] In both figures, the left panel is the distribution weighted by the phase space and the 
Bose-Einstein distribution functions, the middle panel is weighted by the 23 contribution to r\ (denoted as 7723, which 
is the T] analogy of the second term in Eq. (HH)) with the "exact" matrix element, and the right panel is similar to the 
middle one with the GB matrix element. 

We first look at the distribution in the fluid local rest frame in Fig. 3] The left panel plots do not depend on the 
interaction and hence is as independent. The distribution has Q'P&ak ^ ^^^ ^^^ ^^^^^^ variation ae is about the same 
size. In the middle panel, the 7723 weighted distribution with the "exact" matrix element, on the other hand, has 
Qpeak ^ ^/aj at small as, while (B) is significantly bigger and ae is close to its value in the left panel. In the right 
panel, where the GB matrix element is used, Qv^'^^ is still close to be proportional toy^oj at small as, but the angle 
is about twice as big as the "exact" case. 

The distribution in the 23 collision CM frame shown in Fig. [S]has a similar behavior as that in the fluid local rest 
frame but the angles are in general much larger. 

The above analysis suggests that the GB formula, which takes the soft gluon bremsstrahlung limit in the CM frame, 
still has some near collinear splitting behavior in the fluid local rest frame. It is curious what the nature of the long 
tail is. We will leave it for future investigation. 
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Figure 6: Universal curves for (^/s (this work), 77/s |20l | and their ratio. These curves are universal and suitable for a general 
SU{Nc) pure gauge theory. 
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Figure 7: Left panel: ri/s for a pion gas [54] and a gluon plasma with LQCD (6CH62| | and perturbative QCD (20|, together with 
77/s extracted from RHIC elliptical flow (^2): I [13 | and II |13 |. The arrow below the line of "Hydro+«2 data II" indicates that 
it is an upper bound. Right panel: ("/s for a massive [7ll| and massless [sj pion gas (the massive pion curve denoted as "pion 
gas" is expected to merge to the massless pion one when T > m^/3) and a gluon plasma with LQCD [7g | and perturbative 
QCD (this work). The sum rule result is for Nf — 3 [16]. 
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E. More aspects 



In Fig. [6l our results for C/s and rj/s {jj/ s is computed in Ref. [20|) using the "exact" matrix element for the 

23 process are shown in the left panel and their ratio C,/r] in units of a^ and (l/3 — c^) in the right panel. As we 
emphasize in Sec. IIV Al these are universal curves suitable for a general SU{Nc) pure gauge theory. 

The external gluon mass TOoo is included in the entropy density s here, but it is a higher order effect and numerically 
very small at small a^. In the range where perturbation theory is reliable (a^ < 0.1), ^ is always smaller than -q 

by at least three orders of magnitude. One can see that our result of C/ry agrees with 15 (l/3 — c^) of Weinberg 
parametrically [53], and it is rather close to the LL one in Eq. (P5|) . 

In Fig. [71 we have plotted rj/s vs. T /Tc and C,/ s vs. T /Tc for QCD with various number of light quark flavors Nf 
(and different T^s are used in different systems) at zero baryon chemical potential. In rj/ s, T/T^ <C 1, the QCD result 
is calculated by the pion gas system using the Boltzmann equation [5 41 ] (the kaon mass is more than two time bigger 
than Tc — too heavy to be important for T/Tc <C 1; for other calculations in hadronic gases, see |55l - l59j ). The T /Tc > 1 
result is for gluon plasma using lattice QCD (LQCD) |60l - l62| (see [63] for a recent review; for a lattice inspired model 
around T^ see e.g. Ref. [64|). This result has assumed a certain functional form for the spectral function and hence 
has some model dependence. Note that in this temperature region, there might be anomalous shear viscosity arising 
from coherent color fields in the early stage of the QGP [63. We have also shown the value of rj/s extracted from 
the elliptic flow (^2) data of RHIC using hydrodynamics: rj/s = 0.1 ± O.l(theory) ± 0.08(experiment) [13] (denoted as 
"Hydro+W2 data I") and -q/ s < 5 x l/(47r) [I3 (denoted as "Hydro+W2 data 11"). And we have assigned a conservative 
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temperature range T = 0.24 ± 0.10 GeV that covers the initial and final temperatures in the hydrodynamic evolution 
{Tf ^ 0.14 GeV, T, < 0.34 GeV). 

For T/Tc ^ 1, we use the perturbative result of the gluon plasma with the 22 and 23 processes in the Boltzmann 
equation [20] and the standard two-loop renormalization (the scheme dependence is of higher order) for the SU(3) 
pure gauge theory 

where ^o = ll/(167r2) and^i = 102/(167r2)2. Fitting to lattice data at 1.2 < T/Tc < 2 yields n ~ I.Utt, Ajg^ ~ 261 
MeV and T^ ~ 202 MeV [eg. When T/T^ ~ 3.3 and 42, a^ = 0.2 and 0.1, respectively. If rj/s above T^ is dominated 
by the gluon contribution so the gluon plasma result {Nj = 0) is close to that of Nf = 3 QCD, then Fig. [7| shows 
that rj/s might have a local minimum at Tc [54. .67.-.69J. 

In (/s, T/Tc <C 1, the QCD result is calculated by the Boltzmann equation for massless [701 and massive |71| 
(also in Ref. [72l |) pions. For massless pions, C/s is increasing in T since it is expected when the pion self-coupling 
vanishes (or equivalently the pion decay constant /^ — > oo), C, also vanishes. Thus, the dimensionless combination 
(,/ s ex (T I ]t^Y ^ where z is some positive number. For massive pions, however, the 24 process (23 not allowed by parity 
conservation), which is slower than 22, controls the time scale for the system to go back to thermal equilibrium. At 
low enough T, this time scale is very long since there are not many pions energetic enough to collide and produce four 
pions. If the time scale is longer than that of the fire ball expansion at RHIC, the elastic scattering calculation of 
Ref. |55l. P7a. p74| (see also JTal ) might be more relevant phenomenologically. For TjTc ^ 1, lattice QCD calculation of 
a gluon plasma is shown [7a| together with the sum rule result with Nf = 3 [ISl . [la | . Both of them have some model 
dependence on the shape of the spectral function used. This issue was discussed extensively in Refs. |77h79| which 
inspired Ref. |76] to include a delta function contribution to the spectral function which was missed in the earlier 
result of Ref. [80] . The same delta function will modify the sum rule result [iSl [iq as well. This is yet to be worked 
out. 

For T/Tc ^ 1, the perturbative result of the gluon plasma calculated in this work is shown. We see that although 
C/s for the massive pion case is decreasing in T for small T. It should merge to the massless pion result when the 
pion thermal energy ~ 3T is bigger than ttIt^. Thus, it is still possible that C/s has a local maximum at Tc as in some 
model calculations |68ll8lll83] provided there is no much difference between the Nj = and Nj — 3 results above Tc- 

It is very interesting that (> rj for the gluon plasma just above Tc- This suggests a fluid could still be perfect 
without being conformal, like the AdS/CFT model of Ref. [81j. Finally, it is intriguing that rj/s might have a local 
minimum at Tc and ^/s rnight have a local minimum at Tc- However, despite there are many other systems exhibiting 
this behavior for rj/s |58l . 1671 ISSl . |84| , there are counterexamples showing that it is not universal (Sa - ISSJ . 

V. CONCLUSIONS 

We have calculated the shear and bulk viscosity of a weakly interacting gluon plasma with 22 and 23 collisional 
processes and a simple treatment to model the LPM effect. Our results agree with the results of AMY and ADM 
within errors. By studying the 23 contribution to 77, we find that the minimum angle among the final state gluons 9 
has a distribution that is peaked at 9 ^ y^, analogous to the near collinear splitting asserted by AMY and ADM. 
However, the average of 9 is much bigger than its peak value, as its distribution is skewed with a long tail which 
is worth further exploration. The same 9 behavior is also seen if the 23 matrix element is taken to the soft gluon 
bremsstrahlung limit in the CM frame. This suggests that the soft gluon bremsstrahlung in the CM frame still has 
some near collinear behavior in the fluid local rest frame. We also generalize our result to a general SU{Nc) pure 
gauge theory and summarize the current theoretical results for viscosities in QCD. 
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dation of the Shandong province under grant ZR2010AQ008. JD is supported in part by the Innovation Foundation 
of Shandong University under grant 2010GN031. 
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Appendix A: Soft gluon bremsstrahlung 

In this appendix, we give the details of the derivation of the GB formula or the matrix element for the soft gluon 
bremsstrahlung. We work in the CM frame of the initial or final state where the longitudinal direction is defined as 
that of pi or p2. The conditions for the soft gluon bremsstrahlung are: s ^ Pix and kx ^ ygr (or s — )■ oo and ?/ — >■ 0). 
This means the energy of the bremsstrahlung gluon, say E^, is much smaller than other two gluons, E^ <C E^jE^. 

It is convenient to use the Mandelstam-like variables defined as 

S = {Pl +P2f, t = {pi -Pzf, U = (pi -PaY, 
s' = {P3+P4:f, t' = {P2 -Pif, u' = {p2 -Pzf, 

T,,^{h + k,)^, (z = 1,2,3,4). (Al) 

Here we assume all gluons are massless, so we obtain 

(12)-^, (13) = -^, (14) = -|, 

u' , t' ^ s' 

(23) = --, (24) = --, (34) = -, 

(15) = ^, (25) = ^, (35) = ^, (45) = ^. (A2) 

Using light-cone variables in Eq. (jl3[) and taking the limit s — >■ oo or s 3> p|r, we have 



t = 


= -:; (qr-kr)^ 


u s 


a -s, 


1 

s ? 


- (l-y)s. 


t' . 


= -Qt-, 


u' . 


-- -(l-2/)s, 


\5 " 


. 4/y, 


^25 " 


a 2/s, 


_, 


(kr - y^T? 


35 ^ 


{^-y)y 


45 " 


s 2/s. 



(A3) 
We see that t,t' ^T.^ii = 1,2,3,4) are small. In evaluating |Afl2^345|^ we denote (12345) = 1/[(12) (23) (34) (45) (51)], 
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and we can evaluate all quantities in the denominator of Eq. ([5]). 





Sw's'ri5T45 


(12354) = 
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su'uT-ihT^'i ' 


(12435) ^ 


25 


ss't'TisTas' 


(12453) = 
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stt'T:i5TA5 ' 


(12534) ^ 
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ss'uT25n5 ' 


(12543) = 


1 


ss'msTiz ' 


(13245) = 


1 


u'tt'Ti5Ti5 ' 


(13254) = 


1 


Uu'tT25T45 ' 


(13425) = 


1 


s'tt'T,5T25 ' 


(13524) = 


1 


utt'T25T35 ' 


(14235) = 


1 


uu't'Ti^Tss 


(14325) = 


1 


uu's'Ti5T25 ' 



(A4) 

where we have factored out 2^. Note that other permutations which do not appear are given by the identity (12345) — 
(15432). Then we can collect the most singular parts involving t,t',Ti<^(i — 1,2,3,4) in the denominator and obtain 



|Mi2^345|' 



1111 



stt'TssTis u'tt'Ti5T45 s'tt'Ti5T25 utt'T25T35 ' 

1111 



ss't'TisTas ss'tT25T45 uu'tT25T45 uu't'Ti^Ts^ 

1 (i-y)' , 11 y' 



s'^ql 



(qt - kr)'^ (kr - yqr)^ (qt - kr)'^ k^ (ky - yq^rykj 



(A5) 



One can see that the matrix element squared has singularities from three poles at k'^ = 0, (qr ~ k^)^ = 0. For the 
soft limit, kx 3> yqx, this can be realized by setting y — >■ 0, we obtain 



which reproduces the GB formula. 
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Appendix B: Error Estimation 



The error bands of 77 and C shown in Figs. [T]and|3]are based on the estimation of the following errors: 

(a) HTL corrections for the 23 process: In the 22 process, if we replace the HTL scattering amplitude of Eq. ([7]) 
by that of Eq. ([9]) with niD as the regulator, then the 22 collision rate is reduced by ^ 30% for as ~ 0.005-0.1. At 
smaller as, the effect becomes smaller and eventually becomes negligible at as — 10"*. The reduction arises because 
the HTL magnetic screening effect gives a smaller IR cut-off than rriD . Analogously, using rriD as the regulator in the 
23 process tends to under-estimate the 23 collision rate and gives a larger rj and (. 

(b) LPM effect: Our previous calculation on r] using the Gunion-Bertsch formula shows that implementing the 
mo regulator gives a very close result to the LPM effect [42^. Thus, we will estimate the size of the LPM effect by 
increasing the external gluon mass rUg from ttIqc to mo. 
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(c) Higher order effect: The higher order effect is parametrically suppressed by 0{y/al), but the size is unknown. 
Computing this effect requires a treatment beyond the Boltzmann equation |25,] and the inclusion of the 33 and 24 
processes. We just estimate the effect to be ^/al times the leading order which is ~ 10% at as = 0.01. (Note that 
we estimated the higher order effect to be 0{as) suppressed in Ref. [43|. But since the expansion parameter in finite 
temperature field theory is g instead of g^ , we enlarge the error here.) 

Combining the above analyses, we consider errors from (a) to (c). To compute a recommended range of Q (the 
range of r/ is computed analogously), we will work with the R22 and i?23 collision rates defined as 

-^23 = C23, 
{R22 + R23) = C22+23, (Bl) 

where C23 is the bulk viscosity for a collision with the 23 process only. Using HTL instead of m^) for the gluon 
propagator enhances the 22 rate by a factor of 

S^^^^mR. (B2) 

^22(MD) 

We will assume that the same enhancement factor appears in 23 rate as well, such that 

^^mR ^ s. (B3) 

R23{MD) 

On the other hand, the LPM effect is estimated to suppress the 23 rate by a factor of 

^^R23iLPM)_ (34^ 

R23{MD} 

Combining the estimated HTL and LPM corrections to the 23 rate, the 22+23 rate is likely to be in the range 
[R22 + R23,R22 + 7<5i?23], while the higher order effect gives ±y/a7(-R22 +^23) corrections to the rate. Without 
further information, the errors are assumed to be Gaussian and uncorrelated, the total rate is 

R22 + ^^^i?23) ± (^^^^23) ± V^ {R22 +R23), (B5) 

and the recommended upper (C+) and lower (C-) range for C are 

C± = 1^2 (^^) 

R22 + ^R23)t^[^R23) +as{R22+R23f 

The (± values are shown in the left panel of Fig. |3l 
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